Robust quantum memory via quantum control 
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Abstract: A general scheme for building a quantum memory by transferring quantum infor- 
mation to an essentially decoherence-frec memory transition using quantum control is presented 
and illustrated by computer simulations. 
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Quantum computation [1] has been a fruitful area of research lately. Some of the most promising schemes 
involve encoding qubits into ions and neutral atoms in high Q optical cavities. One of the greatest limitations 
to such schemes is the decoherence which occurs at the optical transitions. This decoherence is the limiting 
factor in determining the temporal length of a sequence of pulses to perform a given computation. Many 
schemes have been suggested to overcome this decoherence, especially quantum error correction [2], which 
uses redundant information to compensate for the losses, and decoherence- free substates [3], which use 
combinations of states which are robust against decay because of quantum interference. In this paper we 
employ Lie group decompositions [4] to derive a promising scheme for a quantum memory, with the idea 
being to transfer quantum information from the important channel for quantum information processing to 
a "memory transition" which holds the quantum information. 

Our scheme can be illustrated with a simple example. Consider a four-level atom as depicted in figure 1 
with two degenerate ground states l^i) and I52) and two non-degenerate excited states |ei) and |e2). The 
ground states might be for example Zeeman sublevels of an atom such as Rubidium, although our scheme 
is far more general than this. We assume that the \gi) — |ei) transition is an optical transition that forms 
part of a quantum information processing scheme. After performing some quantum logical operations we 
wish to protect the quantum information stored in this transition in a long-lived system by applying a 
series of Gaussian control pulses (derived from optical fields) to transfer the quantum information onto the 
decoherence-free \gi) — \g2) transition. Although the transition pairs l^i) — |ei), 1^2) — and \gi) — |e2), 
I52) — 1^2) are at equal frequencies, they can be individually addressed through the choice of appropriate 
field polarisations, allowing complete controllability of the system. 

Formally, the problem can be stated as follows. We wish to map the density matrix p representing the 
(initial) state of the system, whose elements are Pgig^, Peiei, Pgiei and P(,2S2 — Pe2e2 ~ Pgig2 — Pgie2 = 
^9261 = Pg2e2 = 0, onto a density matrix p' such that p'g^g^ = Pgigi, P'g2g2 Peiei, p'g^g2 = Pgiei and 
Peiei = Pe2e2 = Pgiei = Pgie2 ~ Pg2e\ = ^9262 = by applying a scquencc of simple control pulses. Because 
the populations and coherence have been mapped to degenerate energy levels with no allowed transitions 
between them, these states will be extremely long lived. The quantum information can then be returned to 
the optical transition simply by using the reverse of the quantum control scheme. In order to realize the 
mapping of p onto p' we find a unitary operator U such that p' = UpW and decompose the operator U into 
a product of simple unitary operators, each realized dynamically by applying a control field 1 (2) (3) that 
drives the \g2) — {62} (Iffi) — |e2)) (l5i) ~ ki)) transition, respectively. Concretely, note that we have 
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Furthermore, we can express J7 as a product U = V1V2V3V2V1 where 
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Fig. 1. Energy level diagram 



Time 

Fig. 2. Control pulses and evolution of the populations and coherence 



Observing that [4] 
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we see that the operator Vi can be dynamically realized by applying a pulse fi{t) = Ai(i)e*"s2=2+'^/2 with 
total pulse area J\Ai{t)\dt — 7r/(2(ig2e2) ^-nd appropriate polarisation, where ujg2e2 is the frequency and 



the absorption oscillator strength of the \g2) 
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1 62) transition. Similarly, 
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shows that V2 [V3] can be dynamically realized by applying a pulse f2{t) ~ A2(t)e*'^si=2+'^/2 [f^{t) ~ 
A3(<)e''^''i=i+^/^] with appropriate polarisation and total pulse area / 1^2(^)1 dt = tt / {2dg^c,^) [/ |A3(i)| dt = 



7r/(2(igjej)] where u> 
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Hence, generation of U requires a sequence of five pulses, where the first pulse drives the transition \g2) — \ e2) 
and has pulse area tt/ {2dg.-^e2)] the second pulse drives the transition \gi) — |e2) and has pulse area tt/ {2dg-^e2)] 
the third pulse drives the transition \gi) — |ei) and has pulse area 7r/(2dgjei); the fourth pulses drives again 
the transition \gi) — \e2) and has pulse area n / (2dg-^e2) ^-^^ the last pulse drives again the transition \g2) — {£2) 
and has pulse area -k I {2dg^e2)- Observe that only the total area of each pulse matters and hence the precise 
pulse envelopes can be adjusted to suit experimental constraints. 

Figure 2 shows the pulse sequence and evolution of a system initially in the superposition state + 
2|ei))/\/5. Note that initially Pg^g^ ~ 0.2, Peiei — 0.8, Pg^ei ~ 0.4 and all other matrix elements of p are 
zero. Observe that at the final time we have indeed p'g-^g-^ = 0.2, p'g^g^ = 0-8, Pgiga = 0.4 and all other 
elements of p' are zero. The time unit in all plots is 1/^2, where is the Rabi frequency of the transition 
l5i>-|ei). 
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